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Abstract 

The  collection  of  high-frequency  (455  kHz)  sides¬ 
can  data  from  tactical  sonars  used  in  mine  counter¬ 
measures  (MCM)  is  typically  for  the  purpose  of  ob¬ 
ject  and  clutter  identification.  The  large  amounts  of 
high-frequency  liigh-resolution  sidescan  imagery  are 
also  a  potential  source  of  data  for  environmental  bat- 
tlespace  characterization.  The  high  resolution  of  the 
455  kHz  sonar  imagery  also  provides  much  informa¬ 
tion  about  the  surficial  bottom  sediments,  however 
their  acoustic  scattering  properties  are  not  well  un¬ 
derstood  at  high  frequencies.  Textural  characteristics 
of  the  high-resolution  imagery  are  therefore  used  to 
represent  the  different  scattering  properties  of  various 
sediment  types.  The  fusion  of  this  data  with  collat¬ 
eral  seafloor  information  from  other  tactical  sensors 
or  databases,  can  be  used  for  classification.  The  op¬ 
erational  goal  is  to  use  the  sidescan-derived  features 
with  a  vertical-incidence  sediment  classifier  to  infer 
sediment  properties  at  angles  off  of  vertical  incidence 
that  can  then  be  used  in  mine  burial  prediction  mod¬ 
els,  acoustic  performance  models,  and  survey  plan¬ 
ning. 

1  Introduction 


Sediment  acoustic-scattering  properties  at  high  fre¬ 
quencies  is  an  area  of  much  interest  and  active  re¬ 
search  [1], [2], [3], [4], [5], [6].  For  tactical  purposes,  the 
probability  of  detection  of  objects  on  the  seafloor 
is  highly  dependent  on  background  scattering  prop¬ 
erties,  particularly  in  cluttered  environments.  The 
high-resolution,  high-frequency  sidescan  sonars  are 
very  well  suited  to  the  task  of  object  detection.  The 
high-resolution  imagery  of  the  acoustic  response  of 
the  seafloor  also  provides  a  large  amount  of  data  of 
opportunity  for  mapping  the  seafloor.  However,  since 
the  high-frequency  acoustic  sediment  physics  are  not 
well  understood,  a  classification  approach  based  on 


textural  segmentation  is  used.  The  fusion  of  this 
data  with  collateral  seafloor  information  from  verti¬ 
cal  incidence,  lower  frequency  multibeam,  or  other 
tactical  sensors  or  databases,  can  be  used  to  classify 
the  low  grazing  angle  data.  The  objective  is  to  use 
the  acoustic  impedance  information  to  label  the  seg¬ 
mented  image  data. 

The  Klein  5500  sidescan  sonar  operates  at  455  kHz 
and  uses  a  dynamically-focused  aperture  to  produce 
five  parallel  beams  per  ping.  The  beams  have  a  con¬ 
stant  resolution  of  either  10  or  20  cm  along  track,  and 
an  across-track  resolution  of  3.75  cm.  The  Klein  5500 
imagery  used  in  this  study  was  collected  by  the  U.S. 
Naval  Oceanographic  Office  off  the  coasts  of  Corpus 
Christi,  TX  and  Panama  City,  FL.  Impedance  data 
was  derived  from  the  UQN-4  vertical  incidence  sonar, 
the  standard  fathometer  for  all  U.S.  Navy  ships.  The 
impedance  data  were  collected  at  different  times  from 
the  sidescan  data. 

The  segmentation  of  the  data  using  several  texture 
algorithms  is  discussed  in  this  paper,  and  some  results 
for  the  Panama  City  data  set  are  presented.  Some 
preliminary  results  are  presented  on  the  fusion  of  the 
impedance  data  with  texture. 

2  Texture  analysis 

Texture  measures  the  patterns  of  acoustic  intensity 
fluctuations, that  can  be  used  to  discriminate  between 
sediments.  Textural  features  are  derived  from  grey- 
level  cooccurrence  matrices  (GLCM)  [7],  frequency- 
domain  measures,  higher-order  moments,  and  fractal 
dimensions  [8], [9]. 

The  GLC'M-derived  features  are  entropy,  homo¬ 
geneity,  contrast,  and  correlation.  The  texture  fea¬ 
tures  are  constructed  from  the  slant-range  corrected 
data  for  several  adjacent  pixel  orientations,  9  = 
0°,  45°,  90°, and  135°  (in  Cartesian  space)  and  inter¬ 
pixel  distances,  d  =  1,2.  The  cooccurrence  probabil- 
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ities  are  averaged  over  the  four  directions  to  stabilize 
the  variance  due  to  the  high  rate  of  towfish  motion, 
and  changes  in  track  heading  as  well.  These  averaged 
features  are  considered  rotationally  invariant. 

The  higher-order  central  moments  of  the  data 
are  used  primarily  to  incorporate  the  relative  mean 
backscatter.  The  mean,  variance,  skewness  and  kur- 
tosis  are  computed  for  each  block  of  pixels  used  in 
the  analysis.  A  mean  backscatter  parameter  is  de¬ 
rived  from  the  data  by  removing  all  angle-varying  or 
range-dependent  gains.  The  resulting  parameter  is 
adjusted  by  an  arbitrary  constant  that  represents  the 
overall  system  gain,  thus  the  mean  backscatter  is  rel¬ 
ative.  The  inclusion  of  the  first  central  moment,  or 
mean,  as  a  feature  in  the  classifier  makes  the  resulting 
segmentation  a  function  of  mean  backscatter  as  well 
as  texture. 

The  fractal  features  are  constructed  for  the  hori¬ 
zontal  and  vertical  directions  following  Linnett  [8],. 
This  approach  derives  the  fractal  dimension  based  on 
the  brightest  reflectors  and  lowest  reflectors,  produc¬ 
ing  an  upper  and  lower  fractal  surface  respectively,  at 
spatial  resolutions  r  =  1,2,  and  4. 

Texture  features  can  be  described  by  the  spatial 
distribution  of  the  two-dimensional  Fourier  trans¬ 
form,  and  the  speed  with  which  the  FFT  can  be  com¬ 
puted,  makes  it  an  appealing  choice  for  near  real  time 
applications.  The  2-D  spatial  spectra  for  an  area  of 
uniform  texture  contain  information  on  the  direction, 
magnitude  and  wavelength  of  spatial  patterns  in  that 
area.  Sand  ridges  are  a  good  example  of  a  directional 
periodic  texture.  If  the  sand  waves  consisted  of  a  sin¬ 
gle  wavelength,  with  crest  and  trough  of  equal  shape 
and  size,  then  the  spatial  spectrum  would  be  an  im¬ 
pulse  at  a  single  point.  As  the  wavelength  varies,  the 
spectrum  broadens  in  the  wave  direction.  If  the  area 
is  homogeneous  with  random  scatterers  then  spec¬ 
trum  will  be  randomly  distributed. 

To  exploit  the  frequency  coherence  seen  in  sand 
ripple  patterns,  a  measure  of  the  squared  coherence 
between  the  horizontal  and  vertical  frequency  com¬ 
ponents  is  a  useful  metric.  The  2D-FFT  is  performed 
over  the  image  using  smaller  overlapping  windows. 
The  data  are  tapered  with  a  Gaussian  window  with 
large  a  to  preserve  spectral  resolution.  The  spectral 
texture  features  are  derived  by  first  identifying  areas 
of  higher  energy  in  the  spectrum.  The  expected  en¬ 
ergy  for  a  2-D  white  noise  field  is  used.  The  spectral 
energy  of  the  image  is  compared  to  the  expected  en¬ 
ergy.  A  binary  image  is  created  by  setting  a  bit  if 
the  energy  exceeds  a  predetermined  threshold.  The 
expected  noise  field  can  be  modelled  in  any  number 
of  ways,  to  satisfy  assumptions  about  the  noise  sta¬ 
tistics. 


A  set  of  {x,  y}  coordinates  is  generated  for  each  bit 
set  in  each  block  of  the  binary  image.  By  comput¬ 
ing  a  least  squares  fit  between  the  {x,  y}  coordinates 
,  the  slope,  /3,  intercept,  a,  and  goodness  of  fit,  R 2 
are  derived.  .  The  statistics  indicate  whether  there 
is  a  strong  linear  trend  (high  R2  )  or  random  scat¬ 
ter,  (very  low  f?2).  The  slope,  offset  and  R2  values 
are  saved  as  features  for  each  block  of  the  image,  al¬ 
though  the  slope  and  offset  are  not  used  for  texture 
parameters. 

The  set  of  features  is  reduced  to  a  few  parameters 
that  correlate  with  some  known  physical  property. 
Presently,  principal  components  are  used  to  reduce 
the  set  of  features  for  classification.  A  more  opti¬ 
mal  approach  is  to  construct  features  based  on  the 
best  separation  of  seafloor  bottom  types  of  interest, 
provided  there  is  sufficient  ground  truth.  Since  the 
data  are  not  segmented  into  classes  to  begin  with, 
and  there  is  little  or  no  ground  truth  to  base  the 
classification  on,  we  have  to  derive  an  estimate  of 
the  statistics  from  the  data.  A  K-means  clustering 
algorithm  is  used  which  reassigns  pixels  to  different 
clusters  while  trying  to  minimize  the  within-cluster 
variance.  The  algorithm  proceeds  until  the  num¬ 
ber  of  clusters  reaches  a  predefined  threshold.  The 
clusters  each  represent  a  bottom  type,  which  is  un¬ 
known,  i.e.  unsupervised  classification.  Classification 
is  performed  using  a  Bayesian  maximum  likelihood 
test  [10].  The  final  classification  produces  a  segmen¬ 
tation  of  the  image  based  on  the  individual  cluster 
statistics.  The  fusion  of  the  texture  data  with  the 
impedance  data  or  some  ground  truth  will  be  used  to 
identify  the  classes. 

3  Segmentation  approach 

All  quantifiable  range-dependent  effects  are  removed 
from  the  sidescan  data  to  make  the  images  represent 
relative  mean  backscatter.  The  backscatter  analysis 
presumes  a  flat  bottom  and  uses  the  sonar  equation 
to  correct  for  systematic  and  random  processes  that 
affect  the  overall  system  gain.  The  system  is  solved 
for  reverberation  using 

RL  =  SL  —  2 TL  +  Sb  +  10  log  A 

where  RL  is  the  area  reverberation  for  a  unit  area  at 
range  r,and  ffy  is  the  scattering  strength  [11].  The 
scattering  strength  is  expressed  using  Lambert’s  Law, 
which  gives  the  scattering  strength  as  a  function  of 
grazing  angle,  9  as 

Sb  =  10  log10  y  +  10  log10  sin2  9 
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where  10  log10  fi  is  a  constant  approximately  equal  to 
the  average  scattering  strength,  independent  of  graz¬ 
ing  angle. 

The  quantity  10  log10  fi  is  estimated  from  the  data 
and  the  result  is  quantized  to  8  bits  for  image  process¬ 
ing..  After  removing  the  angular  dependence  in  the 
data  and  converting  it  to  10  log  //,the  resulting  image 
is  approximately  Gaussian,  whereas  the  original  im¬ 
age  was  closer  to  a  Rayleigh  distribution.  The  more 
normal  distribution  makes  many  of  statistical  algo¬ 
rithms  more  robust. 

The  relative  mean  backscatter,  10  log10  /x, is  com¬ 
pared  between  several  areas  off  Corpus  Christi  known 
to  be  predominantly  sand  and  mud.  The  mean 
backscatter  curves  for  the  sand  and  mud  areas,  i.e. 
ECHO  2  and  CHARLIE  2/DELTA/ECHO  3,  respec¬ 
tively  are  shown  in  figure  1.  Each  line  represents  an 
average  10  log10  /x  over  a  separate  survey  track  and  is 
plotted  with  time  for  comparison.  The  plots  show  a 
clear  separation  between  the  mud  and  sand  profiles, 
with  3-10  dB  of  separation.  This  information  is  use¬ 
ful  for  discriminating  between  the  two  bottom  types 
even  though  a  physical  model  is  not  available. 


Figure  1:  Relative  mean  backscatter  (dB)  vs.  time 
(sec).  Upper  set  of  lines  are  for  ECHO  2,  which  is 
mostly  sand;  the  lower  set  are  for  ECHO  3,  DELTA, 
and  CHARLIE  2,  which  are  mostly  mud. 

The  10  log  (j,  values  are  corrected  for  slant  range 
and  quantized  to  8-bits  for  texture  processing.  The 
GLCM,  moments,  spectral,  gradient,  and  fractal  fea¬ 
tures  are  computed  and  reduced  to  a  set  of  principal 
components.  The  mean  backscatter  is  included  as  a 
feature  in  the  first  moment. 

An  example  of  the  frequency  domain  features  is 
shown  for  the  image  in  figure  2.  The  seafloor  im¬ 


aged  withe  the  Klein  5500  near  Panama  City,  FL  is 
composed  of  sand  ridges  with  an  area  of  shell,  sand, 
and  coral  debris  at  the  bottom  of  the  image.  The 
2D-FFT  is  applied  to  each  block  within  the  image, 
and  results  are  shown  in  figure  3.  The  vertical  and 
horizontal  components  of  the  sand  ridges  are  corre¬ 
lated,  producing  the  linear  pattern,  and  thus  a  high 
r 2  value  when  a  linear  fit  is  applied.  The  results  are 
shown  in  figure  4.  The  image  on  the  left  shows  the 
r 2  values,  and  the  right-hand  image  shows  the  slope 
of  the  line.  The  image  of  the  r2  values  shows  values 
greater  than  0.65  for  the  sand  ridges,  and  much  lower 
values,  less  than  0.20,  for  the  sand,  shell  and  coral  de¬ 
bris.  The  slope  image  shows  the  direction  of  the  sand 
waves  to  be  0.5  to  1.0,  or  just  near  45°. relative  to  the 
horizontal  axis.  The  slope  angle,  corrected  for  ship’s 
heading,  gives  the  direction  of  the  wave  pattern. 


100  200  31  «  500  000 

K5  335-234225  6500:8000 


Figure  2:  Klein-5000  sidescan  image  from  Panama 
City  Survey  Area  BRAVO  (K5-335-234225)  shows 
sand  ridges  at  top,  nadir  at  left. 

Texture  features  were  computed  for  the  two  bot¬ 
tom  types,  sand  ridges  and  an  area  of  shells,  sand 
and  coral  debris,  seen  in  figure  2.  In  figure  5  a  scat¬ 
ter  plot  is  shown  of  contrast  vs.  entropy  for  all  the 
points  in  each  bottom  type.  The  plot  shows  there 
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Figure  3:  Results  of  2-D  (64x64)  FFT  over  adjacent 
image  blocks.  Linear  features  indicate  presence  and 
orientation  of  sand  ridges. 


Figure  4:  Example  of  spectral  parameters:  (a)  direc¬ 
tional  strength  R2  ;  (b)  linear  slope 


is  separation  as  indicated  by  the  line,  although  some 
overlap  is  present.  The  scatter  plot  of  three  texture 
features,  entropy,  contrast,  and  r2,  for  each  bottom 
type  is  shown  in  figure  6.  The  separation  between 
these  two  bottom  types  based  on  the  three  features 
gives  a  better  separation,  as  expected.  A  test  sta¬ 
tistic  can  be  constructed  as  a  linear  combination  of 
the  three  or  more  texture  features  and  used  to  clas¬ 
sify  the  pixels  into  one  of  the  two  classes.  Principal 
components  are  used  to  form  linear  combinations  of 
the  22  texture  measures  to  two  or  three  uncorrelated 
features,  which  are  then  used  in  classification. 
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Figure  5:  Scatter  plot  of  texture  features  entropy 
and  contrast  for  areas  of  sand  ridges  and  and 
sand/shells/coral  debris.  Line  indicates  where  a  dis¬ 
criminant  function  might  be  used  to  classify  based  on 
these  two  features  alone.  Overlap  indicates  existence 
texturally  similar  subareas. 

In  practice,  the  bottom  types  are  not  known  a  pri¬ 
ori  so  a  statistical  test  between  blocks  of  data  is  used 
to  build  homogeneous  statistical  classes.  These  statis¬ 
tics  are  used  as  a  basis  for  segmentation.  Segmenta¬ 
tion  is  performed  using  either  a  maximum  likelihood 
statistic  or  a  Poisson-scattering  based  rule  [12].  The 
classification  result  for  the  image  in  figure  2  using  the 
maximum  likelihood  method  with  the  first  two  prin¬ 
cipal  components  is  shown  in  figure  7.  Results  show 
that  near-nadir  results  are  not  as  good  for  the  sand 
ridge  discrimination  because  of  the  low  contrast  that 
results  from  the  higher  grazing  angle.  However,  some 
features  are  more  robust  to  contrast  than  others,  sug¬ 
gesting  that  a  different  feature  set  should  be  used. 
Also,  the  use  of  training  data  would  be  useful  for  op- 
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timum  feature  selection.  For  most  of  the  off-nadir 
data,  these  results  and  others  [13],  [6],  [12]  ,[1] ,  [14] , 
show  that  several  seafloor  bottom  types  can  be  dis¬ 
tinguished  based  on  textural  characteristics. 
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Figure  6:  Scatter  plot  of  texture  features  (en¬ 
tropy,  contrast,  r2)  from  images  of  sand  ridges  and 
shell/sand/coral  debris. 
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Figure  7:  Texture  segmentation  of  Panama  City 
sidescan  image  shown  in  Figure  1.  Gray  indicates 
sand  ridges,  white  indicates  shell,  sand  and  coral  de¬ 
bris. 

Figure  8  shows  an  overlay  of  impedance  values  to 
a  texture  classification  of  some  of  the  Panama  City 
data.  There  are  two  texture  classes  represented  as 
gray  and  white.  While  the  gray  areas  correlate  with 
higher  impedances  (>  3.0),  the  white  areas  corre¬ 
late  with  lower  impedances(<  3.0)  in  some  areas,  and 
higher  impedances  in  other  areas.  The  gray  area  with 
high  impedances  is  indicative  of  coarse  sand  in  this 
area.  The  white  area  with  lower  impedance  correlates 


with  finer  sand,  while  the  white  area  (upper  left)  with 
higher  impedances  are  likely  areas  of  coarse  sand, 
with  some  coral  debris  and  perhaps  gravel.  There 
is  also  the  possibility  that  the  texture  is  identifying 
a  thin  layer  of  surficial  sediments  over  a  layer  with 
different  sediment  properties  that  can  be  detected  by 
the  low-frequency  sonar,  but  not  the  high-frequency 
sidescan. 

4  Conclusions 

Segmentation  of  sidescan  data  is  useful  for  seafloor 
classification  for  several  reasons.  The  sidescan  data 
are  used  by  the  Naval  Oceanographic  Office  to  up¬ 
date  their  seafloor  bottom-type  databases,  a  manual 
task  best  performed  by  experienced  seafloor  analysts. 
The  use  of  the  texture-segmentation  results  in  a  data 
fusion  approach  to  seafloor  classification  using  col¬ 
lateral  data,  can  be  used  to  infer  sediment  properties 
useful  for  tactical  purposes.  Presently  the  Klein  sides¬ 
can  data  are  being  combined  with  the  UQN-4  vertical 
incidence  data,  from  which  acoustic  impedances  are 
estimated.  The  texture-segmented  data  will  be  clas¬ 
sified  with  acoustic  impedances  from  overlying  tracks 
of  UQN-4  data,  providing  wider  bottom  coverage. 

Several  textural  features  applied  to  high-resolution 
high-frequency  sidescan  data  have  been  presented. 
The  use  of  a  combination  of  features  has  shown  that 
good  discrimination  is  possible  for  areas  of  well  de¬ 
fined  sand  ridges  and  other  textures.  Segmentation 
results  need  to  be  evaluated  with  the  impedance  data. 
In  cases  where  a  texture  shows  different  impedances, 
direct  bottom  sampling  would  be  useful  for  evaluating 
discrepancies. 
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